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I.  INTRODUCTION 

Simulation  is  a  tool  used  in  the  design  ot  new  missile  systems  and  in  the 
modification  or  evaluation  of  existing  systems.  A  missile  simulation  allows 
the  engineer  to  actually  try  out  his  design  without  the  expense  of  building 
and  flying  an  actual  missile.  Many  simplifications  are  made  of  course,  but  a 
great  deal  can  be  learned  through  simulation  with  a  substantial  savings  in 
time  and  expense.  These  simulations  are  usually  put  together  with  emphasis  on 
a  particular  aspect  of  the  design  such  as  guidance,  controlability,  component 
evaluation,  etc.,  and  thus  vary  widely  with  no  single  simulatior  being  a 
general  purpose  analysis  tool  that  is  good  for  every  consideration.  A  central 
part  of  each  missile  simulation  is  the  missile  dynamics  equations  which  govern 
how  the  missile  will  react  to  a  given  force  or  moment.  These  equations  can  be 
implemented  in  a  variety  of  ways  and  are  often  a  source  of  confusion  due  to 
the  differences  and  simplifications  made.  This  report  deals  with  this  issue 
by  developing  several  common  forms  of  the  missile  dynamics  equations.  Methods 
of  implementation  and  some  of  the  possible  simplifications  are  also  discussed, 
and  a  summary  is  included  showing  several  of  the  most  common  forms  and  their 
implementation. 

II.  DEVELOPMENT  OF  THE  EQUATIONS  OF  MOTION 

A.  Coordinate  Systems  and  Newton's  Second  Law 

The  fundamental  law  governing  translational  motion  is  Newton's  second 
law  as  given  below. 

F  ■  d(mv)/dt  (1) 

This  law  is  only  applicable  in  an  inertial  coordinate  frame  and  thus  care  must 
be  taken  to  define  a  suitable  inertial  coordinate  frame.  At  this  point  it 
should  be  noted  that  an  Inertial  frame  can  only  be  defined  by  the  fact  that , it 
is  a  system  in  which  Newton's  laws  hold.  Thus  for  this  article,  an  inertial 
frame  can,  be  taken  to  be  any  frame  such  that  effects  due  to  the  movement  of 
the  frame  are  negligible  for  the  purpose  of  the  simulation.  Two  common  types 
of  coordinate  frames  used  in  missile  simulations  are  coordinate  frames  fixed 
to  the  earth  and  coordinate  frames  that  translate  with  but  do  not  rotate  with 
the  earth.  An  earth  fixed  coordinate  frame  is  a.  coordinate  frame  that  rotates 
and  translates  with  but  does  not  move  with  respect  to  the  earth,  such  as  the 
North-East-Down  system  shown  in  Figure  1. 

One  example  of  the  other  type  is  a  non-rotating  Earth  centered  coordinate 
frame.  This  coordinate  frame  has  its  origin  at  the  center  of  the  earth  and 
translates  with  the  earth  but  does  not  rotate  (see  Figure  2).  That  is  the  . 
origin  of  the  coordinate  system  moves  with  the  earth,  but  the  coordinate  frame 
does  not  rotate  with  respect  to  the  ’’fixed”  stars.  Of  these  two  systems  the 
non-rotating  coordinate  frame  provides  the  "best*  inertial  frame  for  missile 
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simulation  work.  However;  the  earth  fixed  coordinate  frames  can  be  used  in 
certain  cases  such  as  with  non-inertially  guided  missiles  when  the  distances 
traveled  are  small  compared  to  the  earth's  radius  and  the  velocities  are  small 
compared  to  the  velocity  of  escape  from  the  ear“h.  This  will  be  discussed 
more  in  Section  II,  part  D.2. 

Now  assuming  that  we  have  decided  upon  an  inertial  reference  frame, 
Newton's  second  lav  may  be  used  to  find  the  acceleration  of  the  missile  with 
respect  to  this  frame  if  the  forces  acting  on  the  missile  and  the  mass  of  the 
missile  are  known.  The  forces  acting  on  the  missile  can  be  separated  into 
three  main  groups  which  are:  thrust,  aerodynamic,  and  gravitational.  Thrust 
and  aerodynamic  forces  can  be  evaluated  fairly  accurately  with  standard  test 
and  design  procedures  such  as  static  firings  to  obtain  thrust  versus  time 
profiles  for  the  engine  and  wind  tunnel  measurements  to  determine  the 
aerodynamic  characteristics.  Gravitational  forces  can  be  calculated  from  a 
knowledge  of  the  missile's  position  relative  to  the  earth.  The  mass  can  be 
estimated  from  a  knowledge  of  the  missile's  weight  before  and  after  burnout 
(from  test  measurements)  by  using  some  sort  of  relationship  (often  linear)  for 
the  decrease  in  missile  mass  over  the  engine's  burntime.  It  is  not  possible 
to  apply  Newton's  second  law  yet,  however;  as  the  forces  are  only  known  in 
magnitude  and  their  direction  depends  on  the  missile's  orientation  in  the 
inertial  frame.  This  leads  to  the  use.  of  a  missile  body  frame  which  is  a 
frame  fixed  in  the  missile  such  as  shown  in  Figure  3.  The  thrust  and 
aerodynamic  forces  can  be  easily  evaluated  in  both  magnitude  and  direction  in 
this  frame.  All  that  remains  then  is  to  develop  the  relationship  between  the 
Inertial  reference  frame  and  the  missile  body  frame.  To  completely  specify 
the  location  of  the  body  frame  in  the  inertial  reference  frame,  six 
independent  coordinates  are  needed.  These  are  often  referred  to  as  degrees  of 
freedom,  with  a  rigid  body,  in  general,  having  six  degrees  of  freedom  (6-DOF). 
(Note  that  a  rigid  body  is  an  object  that  does  not  bend.)  The  six  degrees  of 
freedom  are  commonly  made  up  of  three  translational  coordinates  and  three 
rotational  coordinates^  The  translational  coordinates  normally  chosen  are  the 
location  of  the  body  frame  origin  in  the  inertial  frame.  Usually  the  missile 
center  of  mass  (c.m.)  is  collocated  with  the  origin  of  the  body  frame.  This 
simplifies  the  calculations  Involved  in  evaluating  the  missile's  acceleration. 
The  rotational  coordinates  must  then  specify  how  the  axes  of  the  body  frame 
are  Inclined  with  respect  to  the  inertial  reference  frame.  This  is  often  done 
using  Euler  angles  which  are  three  angles  that  specify  the  rotational 
orientation  of  the  body  frame;  with  respect  to  the  inertial  reference  frame. 
These  will  be  discussed  in  greater  detail  in  Section  ZI,  part  C. 

B.  Rotational  Equations 

It  is  apparent  from  the  above  discussion  that  to  completely  describe 
the  motion  of  a  missile  in  flight  both  its  translational  and  rotational  motion 
must  be  known.  The  translational  motion  is  governed  by  equation  1  which  is 
the  standard  form  of  Newton's  second  law.  It  would  seem  logical  then  that 
there  would  be  e  similar  equation  governing  rotational  motion,  which  is  In 
fact  the. case.  Thus  the  rotational  equation  of  motion  can  be  stated  as: 

• 

M  -  H 
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(2) 


MISSILE 


Inertial  coordinate 


where  M  is  the  moment  or  torque  acting  on  the  body  and  H  is  the  angular 
momentum.  Note  that  the  reference  point  about  which  the  external  moment  and 
the  angular  momentum  are  calculated  must  be  either  at  the  center  of  mass 
(c.m.)  of  the  body  or  fixed  in  the  inertial  frame.  For  our  purposes  we  will 
assume  that  the  reference  point  (origin  of  the  body  coordinate  frame)  is  at 
the  missile's  center  of  mass.  We  should  also  note  at  this  point  that  the  c.m. 
and  c.g.  (center  of  gravity)  of  the  missile  can  be  taken  to  mean  the  same 
thing  for  our  purposes  due  to  the  relatively  small  size  of  the  missile  when 
compared  to  the  change_in  the  gravitational  field  with  respect  to  distance. 

The  angular  momentum  ”H“  is  given  by: 

H  *  Iw  (3) 

where  I  is  the  inertia  tensor  (3x3  matrix)  and  w  is  the  angular  velocity  of 
the  body.  Note  that  the  elements  of  I  are  not  usually  constant  except  when 
given  in  a  body  fixed  coordinate  system.  Thus  the  uso  of  a  body  Coordinate 
system  simplifies  the  application  of  this  equation.  The  inertia  tensor  can  be 
written  in  matrix  form  as: 


and  the  angular  velocity  vector  can  be  written  in  terms  of  the  missile  body 
rates  (P,  Q,  and  R)  as: 

—  AAA' 

w  -  Pi  +  Qj  +  Rk  (5) 

At  this  point  it  is  important  to  clarify  what  is  meant  by  "referred  to" 
and  "relative  to"  or  “with  respect  to".  The  term  "referred  to"  a  certain 
coordinate  system  means  that  the  vector  is  expressed  in  terms  of  the  unit 
vectors  of  that  system  while  the  terms  "relative  to"  or  "with  respect  to”  a 
particular  system  mean  as  viewed  by  an  observer  fixed  in  that  system  and 
moving  with  it.  Thus  the  same  vector  can  be  referred  to  either  a  moving  or 
stationary  coordinate  system,  but  its  rate  of  change  with  respect  to  time  as 
viewed  by  observers  in  the  two  systems  may  appear  quite  different.  This  is  of 
some  importance  when  the  derivative  of  a  vector  referred  to  a  moving  system  is 
to  be  determined.  Thus  fo t  a  general  vector  "B"  referred  to  a  moving 
coordinate  system  (such  as  a  body  coordinate  frame),  the  time  derivative 
relative  to  an  inertial  coordinate  system  (the  absolute  time  rate  of  change  of 
B)  is:  ' 


B  *  (B)r  +  w  x  B 


(6) 


where  (B)r  is  the  rate  of  change  (with  respect  to  time)  of.  B  as  viewed  by  an 
observer  in  the  rotating  system,  and  w  is  the  angular  velocity  of  the  moving 
frame.  It  is  important  to  note  that  B  is  the  same  vector  in  both  the  inertial 
and  moving  coordinate  frames;  it  is  just  expressed  in  the  unit  vectors  of  the 
moving  system  and  since  these  unit  vectors  could  be  rotating,  the  time  rate  of 
change  of  B  as  seen  by  an  observer  in  the  moving  system  is  not  necessarily  the 
absolute  time  rate  of  change . 


W  *'■•••»  •  m  ft  *  .  *  •-*  *  *  *  '  •  ♦  •  ‘  •  ‘  *  •'#  •  „  *  a  •  ,*JS  ’  i  ■  Jt  <m  »  »  '•  *  «  »  •-  .• 
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Equation  6  can  be  developed  using  the  following  analysis.  Let 


V 

V 
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'B  -  Bjji  +  Byj  +  Bzk 


(7) 


then 


AAA  A  A  A 

B  *  Bjji  +  Byj  +  Bzk  +  Bxi  +  Byj  +  Bzk 


(8) 


AAA' 

where  i,  j,  and  k  are  simply  the  velocities  cf  the  tips  of  the  unit  vectors 
which  are: 


i 


A  —  A 

i  -  w  x  i 


A  _  A 
j-WXj 


(9) 


or 


k  »  w  x  ic 


i 


—  —  A  A  A 

w  x  B  ■  Bxi  +  B  1  +  Bzk 

and  the  derivative  relative  to  the  moving  frame  is: 


(10) 


I 


(B)r  -  Bxi  +  Byj  +  Bzk 


(ID 


The  general  equation  for  the  absolute  rate  of  change  with  respect  to  time  of  a 
vector  referred  to  a  rotating  coordinate  frame  is  then  made  up  of  the  time 
rate  of  change  of  the  vector_relative  to  the  rotating  coordinate  frame  axes 
denoted  by  ()r  plus  an  (w  x  B)  term  to  account  for  the  rotation  of  these  axes 
(unit  vectors)  relative  to  the  inertial  frame. 


Now  returning  to  the  development  of  the  rotational  equations^  equation  6 
can  be  used  to  determine  the  derivative  of  the  angular  -omentum  "H"  referred 
to  the  body  coordinate  frame. 


I 

r 
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H  *  (H)r  +  w  x  H 


where 
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(H)r  -  Iw 


(13) 
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Note  that  equation  8  is  simplified  due  to  the  choice  of  a  body  fixed 
coordinate  frame  (the  derivative  of  the  I  matrix  is  zero  since  it  is  a 
constant  in  a  body  fixed  frame).  Equations  8  and  9  can  be  expanded  as 
follows. 


Xxx 

Xxy 

Xxz 

Iw  * 

Xxy 

lyy 

IyZ 

1 

_  Xxz 

XyZ 

IZZ 

- 

(  XXXR 

+  xxyQ  +  xj 

• 

• 

+ 

(lXyP 

+  IyyQ  +  Ij 

• 

• 

J 

+  dxZP 

+  lyzQ  +  h 

(15) 


(16) 
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and 


w  x  H  -  (HZQ  -  HyR)  i  +  (HjjR  -  HZP)  j  +  (HyP  -  I^Q)  k  (17) 

*  l^xzP  +  xyZQ  +  XZZR>Q  ”  (Xxyp  t  IyyQ  +  IyzR)R]  i 

[(Ixxp  +  *xyQ  +  IxzR)R  *  (Xxzp  +  xyZQ  +  XZZR)P]  j  (18) 


1 

f^xxyP  +  xyyQ  +  IyzR)P  “  (I^x^  +  xXyQ  +  xxzR)Ql  k 

\  Thus 

from  equation  2,  the  moment  acting  on  the  body  is: 

1 

$ 

*  xxxR  +  xxyfQ  "  PR)  *  XXZ(R  +  PQ) 

j 

. 

+  <XZZ  -  Iyy)QR  ♦  IyZ(Q2  -  R2) 

My  -  IyyQ  +  IyZ(R  -  PQ)  +  Ixy(P  +  QR)  (J9 

! 

♦  <xxx  -  IZZ)PR  t  IX2(R2  -  p2) 

*} 

Mz  •  XZZR  +  xX2(P  -  QR)  +  Iyz(Q  +  PR) 

1 

+  (lyy  -  IXX)PQ  ♦  xxy(P2  *  Q2) 

• 

> 

%  '• . 
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*  \ 
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where  1^,  My,  and  Mz  are  the  components  of  the  moment  about  the  x,  y,  and  z 
body  axes  respectively  and  are  governed  by  the  right  hand  rule. 

In  order  to  gain  more  understanding  of  the  possible  simplifications  of 
equations  19,  some  understanding  of  the  inertia  matrix  is  required.  The 
inertia  matrix  can  be  derived  by  considering  the  rigid  body  to  be  a  system  of 
particles  that  do  not  move  with  respect  to  one  another.  The  angular  momentum 
of  the  rigid  body  will  then  be  the  sum  of  the  angular  momentum  vectors  for 
each  of  the  particles.  The  angular  momentum  relative  to  the  point  0  is  then: 


d0  -  l  (Rt  x  MiRi)  (20) 

i 

where  R^  is  the  position  vector  of  the  particle  i  (mass  M^)  relative  to  the 
reference  point  0.  If  0  is  fixed  in  the  body,  then  the  magnitudes  of  all  the 
R^'s  are  constant  and  thus: 

R^  ■  w  x  R^  (21) 

For  a  rigid  body  the  mass  distribution  is  constant,  and  thus  the 
summation  can  be  replaced  by  a  volume  integral  with  mass  pdV  where  dV  is  a 
small  volume  element  and  p  is  the  mass  density  at  that  point.  The  angular 
momentum  Ho  for  a  rigid  body  is  then  given  by: 


Ho  -  /  pR  x  (w  x  R)dV  (22) 

v 


Now  if  the  reference  point  0  is  considered  to  be  at  tke  origin  of  a  cartesian 
system,  with  the  volume  element  dV  located  at  (x,y,z),  the  position  vector  is 
given  by:  ’ 

R  -  Xi  +  Yj  +  Zk  (23) 
and  w  is  given  by 


w  -  Wxi  +Wy  j  +  Wzk 


R  x  (w  x  R) 


* 


A 

k 


[(ZWy  -  YV2)  (XWZ  -  zwx)  (YWX  *  XWy) J 


•  {(Y2  +  Z2)WX  -  XYWy  -  XZWz]  i 

+  [-XYWx  +  (X2  +  Z2)Wy  -  ZYWz]  J 
+  [-XZWx  -  YZWy  +  (X2  +  Y2)WZJ  k 


(24) 


(25) 


(26) 
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The  moments  of  inertia  are  then  defined  as: 


*xx  “ 

/vp(Y2  +  Z2)dV 

hy  = 

/  p(X2  '+  Z2)dV 

V 

(27) 

*zz  “ 

4p(X2  +  Y2)dV 

and  the  products  of  inertia  are  defined  to  be: 

Ixy  -  Iyx  “  -  fv  pXYdV 

Ixz  "I**"  -  /vPXZdv  (28) 

lyz  »  IZy  »  -  fv  pYZdV 

The  angular  momentum  equation  for  a  rigid  body  can  then  be  written  as: 

H  “  (*xxwx  +  *xy^y  +  ^xz^z^* 

+  ( IXyWx  ^  lyywy  ^  ^yi^z^J  (29) 

+  ( *xzwx  +  *yzwy  +  ^zzwz^ 

From  equations  27  it  can  be  seen  that  the  moments  of  inertia  are  always 
positive  and  are  actually  just  a  second  moment  of  the  mass  distribution  with 
respect  to  a  cartesian  axis  system.  The  products  of  inertia  can  be  of  either 
sign  and  can  be  zeroed  by  proper  choice  of  the  axis  system.  In  general  if  a 
three  dimensional  body  has  a  plane  of  symmetry  such  that  the  mass  distribution 
is  a  mirror  image  of  that  on  the  other  side,  then  the  products  of  inertia 
involving  an  axis  perpendicular  to  the  plane  of  symmetry  are  zero  if  the  other 
two  axes  lie  in  the  plane  of  symmetry.  An  example  of  this  is  shown  in  Figure 
5  where  the  missile  shown  is  assumed  to  be  of  constant  density. and  has  two 
planes  of  symmetry  such  that  the  product  of  inertia  terms  are  zero.  Note  that 
missiles  are  not  in  general  of  constant  density,  thus  the  mass  distribution 
must  be  a  mirror  image  on  each  side  of  the  plane  of  symmetry  for  the  product 
of  inertia  terras  to  be  exactly  zero.  Assuming  that  the  product  of  Inertia 

- terms  are  zero  is  a  good  assumption  in  many  cases  though,  for  missiles  having 

two  planes  of  symmetry  such  as  shown  in  Figure  ,5. 

Now  returning  to  equations  19,  it  can  be  seen  that  considerable 
simplification  results  if  the  body  frame  is  chosen  such  that  the  products  of 
inertia  are  zero.  It  was  noted  above  that  a  set  of  axes  could  always  be  found 
such  that  the  products  of  inertia  are  all  zero  and  thus  the  inertia  matrix  is 
a  diagonal  matrix.  The  three  mutually  orthogonal  coordinate  axes  are  known  as 
the  principal  axes  for  this  Case,  and  the  moments  of  inertia  are  known  as  the 
principal  moments  of  inertia.  From  this  it  can  be  seen  that  the  body 
coordinate  frame  for  the  missile  in  Figure  5  is  the  principal  axes  system. 
Equations  19  then  become: 


11 


(30) 


MX  *  *XXP  +  (*ZZ  “  Iyy)QR 

My  ■  IyyQ  +  ( I XX  _  ^ZZ^PP 

Mz  “  *zz^  +  (*yy  ~  *xx^pQ 

which  are  known  as  Euler's  equations  of  motion.  These  are  relatively  simple 
and  are  widely  used  in  solving  the  rotational  motion  of  a  missile.  Note 
that  as  in  the  case  of  the  external  forces,  the  external  moments  are  made  up 
of  aerodynamic  and  thrust  related  moments.  These  can  also  be  estimated  from 
standard  design  and  test  procedures. 

The  rotational  rates  of  the  missile  (P,  Q,  &  R)  can  be  evaluated  using 
equations  30  or  equations  19  as  appropriate.  This,  however;  does  not  provide 
us  with  the  angles  relating  the  body  coordinate  frame  to  the  inertial 
reference  frame  in  a  direct  manner  as  the  angular  rates  (P,  Q,  &  R)  are  not 

the  time  derivatives  of  any  angles  which  specify  the  orientation  of  the  body. 

This  problem  can  be  solved  by  the  use  of  EuLer  angles  or  quaternions,  which 
use  the  angular  rates  to  obtain  coordinate  transformations  from  the  body 
coordinate  frame  to  the  inertial  reference  frame. 

C.  Euler  Angles 

Euler  angles  can  be  used  to  define  the  angular' orientation  of  a 
missile  to  an  inertial  reference  frame.  An  Euler  angle  set  consists  of  three 
angles  and  a  specified  sequence  of  rotation.  In  other  words  to  arrive  at  any 
given  angular  orientation  of  the  missile,  the  missile  axes  can  be  assumed  to 
be  Initially  aligned  with  Che  inertial  reference  frame  axes;  and  then  in  a 
prescribed  sequence,  the  missile  is  rotated  through  each  Euler  angle  about  a 
corresponding  body  axis.  It  should  be  noted  that  there  are  many  different 
Euler  angle  sets  in  use  and  care  should  be  taken  to  define  which  set  is  being 
used  as  different  rotation  sequences  will  usually  give  different  results.  An 
example  of  a  common  set  for  aerospace  engineering  is  shown  in  Figures  6 
through  8.  This  set  consists  of  an  initial  rotation  of  i>  degrees  about  the 
missile  z-.body  axis,  followed  by  a  rotation  of  9  degrees  about  the  missile 
y-body  axis  and  completed  by  a  rotation  of  $  degrees,  about  Che  missile 
x-body  axis.  This  Euler  angle  set  is  also  shown  in  Figures  9 ' through  11  in  a 
pseudo  three  dimensional  view., 


Figure  6'.  Rotation  about  the  z-body  axis  (heading  angle,  ip) . 
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Figure  7.  Rotation  about  the  y-body  axis  (attitude  angle,  9). 


Figure  8.  Rotation  about  the  x-body  axis  (bank  angle.  $). 

In  Figures  6-8  and  9-11,  the  body  axes  are  aligned  with  Intermediate  axes 
systems  denoted  (')  and  (")  after'  the  ^  and  9  rotations  respectively. 

These  are  labeled  so  as  to  facilitate  the  reader's  understanding  but  are  not 
used  otherwise.  Euler  angles  can  be  used  to  specify  any  desired  angular 
orientation  of  the  missile.  Remember  though  that  the  sequence  is  Important 
and  for  the  case  given  another  sequence  such  as  a  9-^ -$  sequence  would 
in  general  give  quite  different  results. 

If  the  inertial  frame  shown  in  Figures  9-11  is  an  earth  fixed  coordinate 
frame  such  as  an  N-E-D  (North-East-Down)  coordinate  frame,  the  Euler  angles 
are  sometimes  referred  to  as  the  heading  angle  (t) ,  the  attitude  angle 
(9),  and  the  bank  angle  (4).  Thus  letting  the  N-E-D  axes  correspond  to 
the  2^  axes  respectively,  the  angle  ^  would, simply  be  the  heading  of 

the  missile  relative  to  north,  the  angle  9  would  be  the  attitude  of  the 
missile  relative  to  the.  north-east  plane,  and  the  angle  a  would  be  the 
lateral  inclination  or  bank  of  the  missile.  Note  that  the  ^-9-$ 
sequence  used  here  could  also  be  referred  to  as  a  yaw,  pitch  and  roll  sequence 
in  missile  terminology. 

The  Euler  angle  rotations  shown  in  Figures  6-8  or  9—1 1  can  now  be  used  to 
develop  a  coordinate  transformation  matrix  that  will  transform  vectors  from 
the  inertial  frame  to  the  body  frame.  Referring  to  Figure  6,  the  coordinate 
transformation  from  the  inertial  to  the  intermediate  ('.)  system  can  be  seen  to 
be : 
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(31) 


~x'~ 

COS  l|/ 

sin 

0 

y* 

m 

-sin  ij) 

cos 

0 

z' 

0 

— 

0 

1 

Note  that  in  the  preceding  transformation,  the  rotation  i3  about  the  Z±  axis 
and  thus  z'  is  collocated  with  Z^.  In  a  similar  manner,  the  rotation  from  the 
intermediate  (’)  system  to  the  second  intermediate  (")  system  can  be  seen  from 
Figure  7  to  be: 


Y"  axis 


Figure  10.  Rotation  about  the  intermediate  Y* . 


Figure  11.  Rotation  about  the  x-body  axis. 
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where  the  rotation  was  about  the  y'  axis  and  thus  the  y'  and  y" 
same.  Then  from  Figure  8,  the  final  transformation  is: 


(32) 


axes  are  the 


[:  - 


0 

cos  $ 
sin  $ 


sin 

cos 


A  [?] 


(33) 


where  the  rotation  wai.  about  the  x"  axis  which  was  also  the  x-body  axis. 
These  transformations  can  be  combined  into  one  transformation  as  follows, 
equation  33  substitute  equation  32  for  (x",y",z"),  then 
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0  -s  9 
1  0 
0  c  9 


x' 

r 

z' 


— s  9 
8  $  *c  9 
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(34) 


(35) 


where  the  sine  and  cosine  terms  have  been  written  in  a  simplified  notation. 
Now  substitute  equation  31  for  (*' ,y* ,*') 


thus 


Phi 

'c  9 

0  -s  9  1 

f  c’4> 

s  i|> 

°1 

.  ;• 

- 

S  <fr  *3  9  c 

$  9  *  *c  9 

Clj) 

0 

Yt  (36) 

* 

UbJ 

.c  $  *s  9  -8 

6  c  ♦  *c  9  J 

L  o 

0 

lJ  1 

LzJ 

» 

'■  1 

V 

m 

c  9  *c 

c  9  *s  \|>. 

-s 

9 

*i 

1 

■  i 

e 

i  »  • 

Yb 

s  <>  *s  9  *c 

s  $  *s  9  *s£ 

s  $ 

*C  9 

Yi 

| 

V 

m 

•  •  c  ♦  *s  t 

+  C  (*  *c4> 

(37) 

i 

i 

Zb 

, 

c  ♦  *s  9  *c 

c  o  *s  9  *s 

c  $ 

*c  9 

Zi 

I 

.  +  s  9  *s  iji 

-  s  *  *c  ,4> 

m 

*  - 

>n  37 

can 

now  be  used  to 

transform  vectors  expressed 

in 

the 

inertial 

'i 

frame  to  the  body  frame.  In  other  words  vectors  referred  ,,to  the  inertial 
frame  (expressed  in  the  inertial  frame  unit  vectors)  are  transformed  so  that 
they  are  referred  to  the  body  frame  (expressed  in  the  body  frame  unit 
vectors).  A  transformation  from  the  body  frame  to  the  inertial  frame  can  also 
be  developed  from  equation  37  by  noting  that  a  transformation  matrix  is  an 
orthogonal  matrix  and  thus  its  Inverse  is  just  its  transpose.  The  Inverse 
transformation  can  then  be  obtained  by  taking  the  transpose  of  the 
transformation  matrix  in  equation  37.  This  can  be  better  understood  by  first 
expressing  equation  37  In  matrix  form  as: 
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*  •* 


[B]  -  [THI] 


(38) 


where  [B]  is  the  body  vector  matrix  ,  [I]  is  the  inertial  frame 

vector  matrix  (X^.Y^.Z^)  and  [T]  is  the  transformation  matrix.  Now  multiply 
by  the  inverse  of  the  transformation  matrix. 


thus 


*1 

*i 

Zi 


t  t 

[T]  [B]  -  [T]  [Tj[I]  -  [I] 


(39) 


[I]  -  [T] 
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c  9  *c 
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c  0  *s  4> 
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+  c  ♦  *c  4> 

-s  a 

s  '  $  *c  8 
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K| 

+  s  $  *s  4> 

c  ♦  *S  0  *8  ||| 

Yb 

-  s  $  *c  4* 

c  $  *c  9 

(40) 


(41) 


Equation  41  can  be  used  to  transform  vectors  referred  to  the  body  system  (such 
as  missile  acceleration,  velocity  or  angular  rate)  to  vectors  referred  to  the 
inertial  system. 


Now  that  it  has  been  shown  how  to  obtain  the  transformation  matrix  when 
given  the  Euler  angles,  the  next  thing  is  to  determine  the  Euler  angles  during 
the  missile  flight.  This  can  be  done  by  determining  the  missile's  initial 
Euler  angles  at  launch  relative  to  the  inertial  frame  and  then  updating  these 
over  the  missile  flight  using  the  angular  rates  of  the  missile. 

Initial  Euler  angles  (before  missile  launch)  can  be  determined  using  the 
definitions  of  the  Euler  angles.  Thus  the  heading  angle  tj>  is  the  angle 
between  the  missile  pointing  direction  and  the  inertial  X  axis  on  the  X^-Y^ 
plane,  the  attitude  angle  6  is  the  angle  between  the  X^-Y^  plane  and  the 
Xj,  axis,  and  the  bank  angle  $  is  the  rotation  of  the  missile  about  the  Xj, 
axis  with  respect  to  the  Xj-Y^  plane.  These  angles  can  be  used  in  the 
transformation. matrix  of  equation  41  to  provide  the  initial  transformation 
from  body  coordinates  to  Inertial  coordinates.  Once  the  missile  Is  launched, 
the  apgular  rates  of  the  missile  (P,  Q,  &  R)  can  be  determined  using  equation 
19  or  30  and  then  related  to  the  Euler  angle  rates  (y;  0,  a)  using: 

P  •  $  -  '  4*  *sin  9 

Q  •  9  *cos  $  +  4>  *cos  9  *sin  $  (42) 

•  • 

R  •  4>  *cos  9  *cos  <>  •*  0  *sin  $ 

Equations  42  were  developed  from  inspection  of  Figure  12,  which  shows  the 
relative  orientation  of  P,  Q,  and  R  along  with  \p ,  0,  and  £•  Solving 
equations  42  for  41 ,  9,  and  i  then  yields: 


Figure  !.?.  Relationship  between  missile  body  rates  and  euler  angle  rates. 


(43) 


'i'  »  vQ*sin  4>  +  R*cos  <J>  )/cos  9 

9  *  Q*cos  0  -  R*sin  <t 

<fc  *  P  +  *sin  0 

*  P  +  (Q*sin  4>  +  R*cos  <fr)*tan  9 


Equations  43  can  be  integrated  to  obtain  and  $  which  can  then  be 

used  to  obtain  the  transformation  matrix  in  equations  41  and  37.  Note  that 
the  equation  has  a  singularity  at  0  ■  90  or  -90. 

0.  Translational  Equations 

1*  Translational  Equations  in  an  Inertial  Frame 

Since  we  now  have  a  way  to  relate  vectors  expressed  in  the  body 
frame  to  the  inertial  frame,  the  translational  equations  of  motion  can  be 
developed.  Equation  1  can  be  rewritten  as: 


The  term  vm  is  due  to  mass  ejection  from  the  rocket  motor  and  is  Included  in 
the  thrust  force.  Thus  equation  45  can  be  written  as: 


F'  ■  ma  (46) 


where  7'  is  the  total  force  on  the  missile  and  is  the  sum  of  the  thrust, 
aerodynamic,  and  gravitational  forces,  a  is  the  acceleration  of  the  missile 
relative  to  the  inertial  frame  and  m  is  the  instantaneous  mass  of  the,  missile. 
The  acceleration  of  the  missile  referred  to  the  body  frame  is: 

Abx  *  FV» 

Aj,y  •  F'y/m  (47) 

Abz  -  F’z/m 

Note  that  the  components  of  acceleration  in  equation  47  are  referred  to  the 
body  frame,  as  F*  is  more  easily  evaluated  in  this  frame.  In  other  words  AjjX, 
A^y,  and  Abz  are  referred  to  the  body  frame  since  F**,.  F*y,  and  ?'z  are 
referred  to  the  body  frame.  '  If  an  accelerometer  triad  was  collocated 
with  the  missile  axis  system,  these  components  would  be  equivalent  to  the 
accelerometer  readings  plus  gravity.  The'  velocity  of  the  missile  can  then  be 
determined  using  numerical  Integration,  but  care  must  be  taken  as  Integrating 
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equations  47  directly  will  not  give  the  correct  answer  since  the  acceleration 
vector  is  referred  to  a  rotating  coordinate  frame.  In  order  to  get  the 
correct  expression  for  the  velocity  derivative,  equation  6  for  the  derivative 
of  a  vector  in  a  rotating  frame  must  be  used*  This  then  results  in: 


Ab  *  (Vr  +  w  *  vb  (48) 

The  angular  rate  (w)  and  the  velocity  (Vb)  of  the  missile  are  referred  to  the 
body  system  and  are  given  by: 

(49) 

(50) 

A  A  A 

where  ib,  jb>  kb  are  the  unit  vectors  in  the  body  frame.  Thus  from  equations 
49  and  50 

w  x  Vb  -  (QW  -  RV'ijj  +  (RU  -  PW)jb  +  (PV  -  QU)kb  (51)  • 

and  equation  48  can  be  expanded  and  combined  with  equation  47  to  give: 

vbx  -  *  V°  “  (QW  ~  HV) 


where  Che  F/rn  terms  are  Che  accelerations  Chat  would  be  Indicated  by  an 
accelerometer  triad  collocated  with  the  missile  body  axis  system  and  “G"  is 
the  acceleration  due  to  gravity  (mass  attraction).  Equation  53  can  be 
integrated  to  obtain  the  velocity  of  the  missile  referred  to  the  body  frame. 
The  coordinate  transformation  of  equation  41  can  then  be  used  to  determine  the 
velocity  of  the  missile  referred  to  the  inertial  frame,  and  as  the  inertial 
frame  is  not  rotating,  these  velocity  comj- :>nents  can  be  integrated  directly 
to  obtain  the  displacements  referred  to  the  inertial  frame. 

Another  method  which  could  be  used  to  determine  the  velocity  and 
displacement  of  the  missile  is  to  use  equation  41  to  transform  the 
accelerations  in  equation  47  to  the  inertial  frame.  The  accelerations  would 
then  be  referred  to  a  non-rotating  frame  and  could  be  integrated  to  obtain  the 
velocity  and  displacement  without  the  addition  of  any  cross  product  terms. 

For  missile  simulation,  the  velocity  of  the  missile  referred  to  the  body  frame 
is  usually  heeded  for  aerodynamic  calculations.  This  can  be  determined  by 
using  equation  37  to  transform  the  velocity  from  the  inertial  frame  to  the 
body  frame. 

It  is  important  to  remember  that  when  we  speak  of  transforming  a  vector, 
we  mean  that  the  vector  is  expressed  using  a  new  set  of  unit  vectors,  and  the 
vector  itself  has  not  been  changed.  Thus  transforming  a  vector  from  the  body 
frame  to  the  inertial  frame  is  just  using  the  components  of  the  vector 
referred  to  the  body  frame  and  the  transformation  equation  (equation  41) 
to  find  the  components  of  the  same  vector  referred  to  the  inertial  frame. 

2.  Translational  Equations  for  a  Moving  Frame 

In  many  cases  it  is  desired  to  obtain  the  location  and  velocity  of 
the  missile  with  respect  to  a  given  point  on  the  earth.  This  is  especially 
important  for  inertially  guided  missiles,  which  need  to  know  their  location 
relative  to  the  earth  in  order  to  arrive  at  the  target  point.-  For  this  case 
an  earth  fixed  system  is  very  convenient  to  work  with,  but.  as  it  is  not  an 
inertial  system  its  movement  must  be  considered.  This  can  be  done  in  the 
development  of  the  missile  equations  of  motion  and  is  shown  in  the  following. 

In  order  to  obtain  the  equations  of  motion  for  a  missile  in  an  earth  ___ 
fixed  coordinate  frame,  equation  46  can  be  written  in  terms  of  the  vector  R^ 
which  gives  the  location  of  the  missile  with  respect  to  the  inertial  frame. 
Thus 


F’  •  mRi  (54) 

where 

Rj.  -  R  +  r  (55) 
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and 


R^  =*  R  +  r 


(56) 


These  vectors  are  shown  in  Figure  13.  Since  _r  is  referred  to  the  rotating 
earth  fixed  frame,  the  second  derivative  of  r  in  equation  56 
can  be  evaluated  using  equation  6. 


r  -  (r)r  +  w  x  r  (57) 

and 

•  •  •  •  •  • 

r  ■  (r)r  +  w  x  r  +  2w  x  (r)r  +  w  x  (w  x  r)  (58) 

Now  using  equation  56  we  h**.ve: 

*  *  *  *  t  • 

Ri  -  (r)r  +  w  x  r  +  2w  x  (r)r  +  w  x  (w  x  r)  +  R  (59) 


where  w  is  the  angular  velocity,  of  the  earth  and  R  is  the  acceleration  of  the 
earth  fixed  coordinate  frame's  origin.  The  equations  of  motion  for  an  earth 
fixed  coordinate  frame  can  now  be  written  as: 


a  »  F*/ra  i-  wx(wxr)  -2wxv-wxr-R  (60) 

where  v  is  the  velocity  and  a  is  the  acceleration  ojf  the  missile  referred  to 
the  earth  fixed  system.  For  an  earth  fixed  frame,  w  is  constant  and  the 
acceleration  of  the  coordinate  frame's  origin  is  only  due  to  the  earth's 
rotation.  Thus  equation  60  can  be  written  as: 


a  •  F'/o  +  w  x  (w  x  R)  •  »*(»*?)*  2w  x  v  (61) 


This  is  sometimes  written  as: 
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a  *  F/m  +  g  +  w  x  (w  x  R)  -  w  x  (w  x  r)  -  2wxv  (62) 


where  the  forces  have  been  split  into  gravitational  (mass  attraction) 
acceleration  "g"  and  other  external  forces  "F“  (aero  and  thrust).  Foj:  most^ 
missile  simulations  the  distance  flown  is  small  enough  that  the  term  w  x  (w  x 
r)  in  equation  62  is  negligible  and  can  be  neglected  without  loss  of  accuracy. 

The  general  equations  of  motion  referred  to  an  earth  fixed  coordinate 
frame  are  then  given  by: 


^  -  Fx/m  +  Gx  -  (Wy  +  W2)Rx  +  WxWyRy  +  WXWZRZ  -  2(WyVz  -  WzVy) 

2  2, 

Ay  -  Fy/m  +  Gy  -  (Wx  +  Wz)Ry  +  WxWyRx  +  WzWyRz  -  2(WZVX  -  WXVZ)  (63) 

2  2 

Az  -  Fz/m  +  Gz  -  (Wx  +  Wy)Rz  +  WxWzRx  +  WyWzRy  -  2(WxVy  -  WyVx) 

where  the  subscripts  x,y,z  stand  for  the  components  along  the  x,y,z  axes  of 
the  earth  fixed  frame.  The  coordinate  transformation  used  should  now  relate 
the  body  and  earth  surface  frames  thus  allowing  the  missile  body  forces  or 
accelerations  to  be  transformed  to  the  earth  fixed  frame.  Equations  63  can  be 
integrated  twice  to  obtain  velocity  and  displacement  of  the  missile  referred 
to  the  earth  fixed  frame  and  the  velocity  of  the  missile  referred  to  the 
missile  body  frame  can  be  obtained  using  the  inverse  coordinate 
transformation. 

A  common  earth  fixed  coordinate  frame  is  the  North-East-Down  system  shown 
in  Figure  1.  For  this  frame  the  earth  radius  is  in  the  down  or  z  direction 
thus  Rx  and  Ry  are  zero  and  Wy  is  zero  as  the  east  direction  is  perpendicular 
to  the  rotation  vector.  Thus  equation  63  can  be  reduced  to: 


Ax  -  Fx/m  +  Gx  +  WXWZRZ  +  2WzVy 


Ay  -  Fy/m  +  Gy  -  2(WZVX  -  Wxvz)  (64) 

2 

Az  -  Fz/m  ♦  Gz  -  WXRZ  -  2WxVy 

Now  noting  that  the  earth's  rot-tion  rate_is  OU000073  rad/sec  and  the  earth's 
radius  is  approximately  6,371,000  m,  the  w  x  (w  x  R)  (centrifugal 
acceleration)  terms  In  equation  64  can  be  estimated  as  follows. 
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2 

VXW2R2  *  -  W  cos(lat)sin(lat)  *  R 

2 

-  -  (0.000073)  *  cos(lat)sin(lat)  *  6,371,000 
■  -  .034  co8(lat)sin(lac)  m/s2  (65) 


2  2  2 
Wx  Rz  -  W  cos(lat)  *  R 

2 

•  .034  cos(lac)  m/s2  (66) 

Thus  the  centrifugal,  acceleration  terms  are  small  and  may  be  neglected  for 
many  cases  depending  upon_the  desired  accuracy  of  the  simulation.  This  also 
justifies  neglecting  tfie  w  x  (w  x  r)  term  In  equation  62  for  cases  where  the 
range  of  the_mls£lle  “r"  Is  much  less  than  the  earth's  radius  "R".  The 
corlolls  (  2w  x  v  )  terms  may  also  be  estimated  If  the  missile  velocity  Is 
known.  If  ve  assume  that  the  maximum  missile  velocity  Is  Mach  2  or 
approximately  670  m/s,  then  the  corlolls  accelerations  will  be  on  the  order 
of  j 

2WV  -  2*0.000073*670  -  .1  ra/s2  (67) 

Thus  corlolls  acceleration  can  also  bn  neglected  for  many  cases ,  but  becomes 
more  important  for  the  faster  missiles.  Earlier  It  was  stated  that  one  of  the 
requirements,  for  assuming  the  earth  fixed  frame  to  be  an  inertial  frame  was 
that  the  velocity  of  the  missile  must  be  much  smaller  than  the  velocity  of 
escape  from  the  earth.  As  the  escape  velocity  is  approximately  11,200  m/s, 
this  can  be  seen  to  be  a  good  approximation.  From  this  simple  analysis.  It 
can  be  seen  that  equations  63  can  be  simollfied  to: 


Ax  "  V“  +  °x 

Ay  -  Py/m  +  Gy  (68) 

A2  -  P2/m  +  G2 

for  many  cases  involving  short  range,  slower  flying  missiles  which  do  not  use 
Inertial  guidance. 

From  the  above  calculations.  It  might  be  thought  that  the  corlolls  and 
centrifugal  acceleration  can  be  neglected  for  all  flights,  but  it  must  be 
remembered  that  even  small  errors  in  acceleration  can  become  significant  if 
they  are  constant  over  a  long  time.  To. get  distance  these  terms  are 
Integrated  twice  which  gives  a  time  squared  multiplier.  Thus  the  distance 
curves  for  a  constant  error  are  parabolic  as  shown  in  Figures  14  and  15. 

The  error  in  displacement  due  to  the  centrifugal  acceleration  at  a  latitude  of 
30  degrees  (calculated  from  equation  66)  Is  shown  In  Figure  14.  In  Figure 
15  the  error  In  displacement  due  to  corlolls  acceleration  is  shown  for 


Figure  15.  Displacement  error  due  to  neglecting  corlolls  acceleration  terms 


missiles  flying  at  constant  velocities  at  0  degrees  latitude.  Note 

that  these,  are  only  approximations  to  give  the  reader  an  idea  of  the  magnitude 

of  the  errors  involved. 

lit.  SUMMARY 

In  order  to  clarify  some  of  the  preceding,  a  short  summary  of  the  possible 
choices  of  equations  is  presented  below.  Note  that  these  are  not  the  only 
options  available,  but  they  are  some  of  the  most  popular.  For  each  of  these 
choices,  a  body  axis  system  and  an  inertial  or  earth  fixed  coordinate  system 
are  needed  and  should  be  defined  first  along  with  thqir  relative  orientations. 
The  moments  and  forces  (or  accelerations)  along  with  the  mass  and  moments  of 
inertia  are  assumed  to  be  known  in  the  body  system. 

A.  Rotational  Equations 

The  rotational  equations  (equations  19  or  30)  are  the  same  for  each 
choice  and  are  used  to  obtain  the  missile's  angular  rates  (P,Q,R)  from  which 
the  transformation  matrix  relating  the  inertial  and  body  coordinate  systems 
can  be  developed.  Thus 

•  •- 
Ixy(Q  -  PR)  +  IX2(R  +  PQ) 

2  2 

-  Iyy)QR  +  IyZ(Q  -  R  ) 

My  -  IyyQ  +  IyZ(R  ~  PQ)  +  I*y(P  +  QR)  (19) 

2  2 

♦  dxx  ",  >  W*  -  P  > 

Mz  -  IZZR  +  IXZ(P  -  QR)  ♦  Iyz(Q  ♦  PR) 

2  2 

,  +  <Iyy  “  WPQ  ♦  Vp  *  Q  > 

or  for  a  missile  with  a  high  degree  of  symmetry,  Euler's  equations  of  motion 
may  be  used: 

**x  "  Ixxp  +  (I*z  ”  Iyy)QR 

My  -  IyyQ  +  (In  -  Izz>PR  (30) 

Mz  “  IZZR  +  (Iyy  “  IXX)PQ 

There  are  several  alternatives  available  for  solving  the  translational 
equations.  These  are  presented  below. 


Ixx*  + 


+  (I 


zz 
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B.  Translational  Equations  in  the  Inertial  System 

One  of  the  simplest  methods  is  to  transform  the  missile  acceleration 
components  from  the  body  coordinate  frame  to  the  inertial  frame  and  then 
integrate  to  obtain  velocity  and  displacement  referred  to  the  inertial  axes 
system.  The  velocity  referred  to  the  missile  body  frame  can  then  be  obtained 
by  using  the  inverse  transformation.  For  this  case  the  body  accelerations 
to  be  transformed  and  integrated  would  be  those  sensed  by  the  missile's 
accelerometers  plus  gravitational  acceleration  due  to  mass  attraction. 


Ax  *  V®  + 

Ay  «  Fy/m  +  Gy  (69) 

Az  «  Fz/m  +  Gz 

The  terms  in  equations  69  are  referred  to  the  inertial  coordinate  frame  where 
the  "F/rn"  terms  are  equivalent  to  the  accelerations  measured  by  the  missile's 
accelerometers  once  they  have  been  converted  to  the  inertial  frame  and  the  "G“ 
terms  are  gravitational  acceleration  (mass  attraction). 

C. "  Translational  Equations  in  the  Body  System 

Another  simple  method  is  to  use  equations  53  to  calculate  the  time 
derivative  of  velocity  referred  to  the  missile  body  axes  system. 

Vbx  *  Fx/m  +  Gx  -  QW  +  RV 

Vby  -  Fy/m  ♦  Gy  -  RU  +  FW  (53) 

Vbz  -  Fz/m  +  Gz  -  PV  +  QU 

These  can  then  be ' Integrated  to  obtain  the  missile  velocity  referred  to  the 
body  frame.  As  displacement  is  usually  desired  in  the  inertial  or  an  earth 
fixed  frame,  the  velocity  components  obtained  are  generally  transformed  to  the 
inertial  frame  and  integrated  directly  to  obtain  the  displacements  referred  to 
the  inertial  frame  axes  system. 


Note  that  Fx/m,  Fy/m,  and  Fz/m  would  be  the  accelerations  measured  by  an 
accelerometer  triad  collocated  with  the  missile  body  coordinate  frame  and  that 
Gx,  Gy,  and  Gz  would  be  the  acceleration  due  to  gravity  (mass  attraction)  along 
these  axes. 

D.  Translational  Equations  in  an  Earth  Fixed  System 

In  order  to  obtain  the  missile's  motion  with  respect  to  an  earth  fixed 
frame,  equations  63  can  be  used. 
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2  2 

^  -  Fx/ffl  +  Gx  -  (Wy  +  Wz)Rx  +  WxWyRy  +  WXWZRZ  -  2(WyVz  -  tfzVy) 

2  2 

Ay  -  Fy/a  +  Gy  -  ( Wx  +  Wz ) Ry  +  WxWyRx  +  VzWyRz  -  2(WZVX  -  Wxvz>  (63) 
2  2 

Az  -  Fz/m  +  Gz  -  (Wx  +  Wy)Rz  +  WxWzRx  +  WyWzRy  -  2(WxVy  -  WyVx) 


or  for  the  specific  case  where  the  earth  fixed  frame  is  a  North-East-Down 
coordinate  frame,  these  can  be  reduced  to: 

Ax  -  Fx/m  +  Gx  +  «xwzRz  +  2Wzvy 

Ay  -  Fy/m  +  Gy  -  2(WZVX  -  WXVZ)  (64) 

2 

Az  -  Fz/a  +  Gz  -  WxRx  -  2WxVy 

Note  that  all  the  terms  in  equations  63  and  64  are  referred  to.  the  earth  fixed 
system.  Thus  for  this  case  Fx/m,  Fy/m,  and  Fz/m  are  the  accelerations  that 
would  be  obtained  from  the  missile's  accelerometers  after  transforming  them  to 
the  earth  fixed  coordinate  frame.  Gravitational  acceleration  due  to  mass 
attraction  is  included  as  a  separate  term.  These  accelerations  can  be 
integrated  twice  to  obtain  the  velocity  and  the  displacement  of  the  missile 
referred  to  the  earth  fixed  frame. 

5.  Translational  Equations  with  the  Earth  Fixed  System  Assumed  to 

be  an  Inertial  System 

A  somewhat  simplified  case  that  works  for  some  of  the  shorter  range, 
slower,  non-inert tally  guided  missiles  is  to  reduce  equations  63  to 


A*  -  Fx/m  +  Gx 

Ay  -  Fy/m  +  Gy  (68) 

Az  -  Fz/m  +  Gz 

by  neglecting  the  centrifugal  and  coriolls  acceleration  terms.  Note  that 
again  the  terms  in  these  equations  are  referred  to  the  earth  fixed  system. 

Thus  for  this  case  Fx/m,  Fy/m,  and  Fz/m  are  once  again  the  accelerations  that 
would  be  obtained  from  the  missile's  accelerometers  after  transforming  them  to 
the  earth  fixed  coordinate  frame  and  "G"  is  gravitational  acceleration  due  to 
mass  attraction.  This  method  is  equivalent  to  alternative  B  with  the  earth 
fixed  system  taken  as  the  inertial  frame.  Thus  equations  68  can  be  integrated 
twice  to  obtain  the  velocity  and  displacement  of  the  missile  referred  to  the 
e«.rth  fixed  frame. 
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F.  Translational  Equations  In  the  Body  System  with  the  Earth 
Fixed  System  Assumed  to  be  an  Inertial  System 

Since  alternatives  B  and  E  are  equivalent  If  the  earth  fixed  frame 
is  assumed  to  be  an  inertial  frame,  an  alternative  equivalent  to  alternative  ■ 

C  is  also  possible  if  the  same  simplflcatlons  used  in  alternative  E  are 
used.  Thus  the  equations  of  alternative  C  are  again  obtained  where  this 
time  the  earth  fixed  frame  is  assumed  to  be  an  inertial  frame. 

Vbx  -  Fx/»  +  G*  -  QV  +  RV 

vby  -  Fy/m  +  Gy  -  RU  +  PW  (70) 

Vbz  -  Fz/m  +  Gz  -  PV  +  QU 

These  can  then  be  Integrated  to  obtain  the  missile  velocity  referred  to  the 
body  frame.  The  velocity  obtained  can  then  be  transformed  to  the  earth  fixed 
frame  and  integrated  to  obtain  the  displacement  referred  to  the  earth  fixed 
frame.  Note  that  Fx/m,  Fy/a,  and  Fz/m  are  the  accelerations  that  would  be 
measured  by  an  accelerometer  triad  collocated  with  the  missile  axis  system  and 
"G"  is  the  acceleration  due  to  gravity  (mass  attraction).  The  terms  in 
equations  70  are  referred  to  the  missile  body  coordinate  frame. 

IV.  CONCLUDING  REMARKS 

The  missile  dynamics  equations  have  been  developed  using  the  basic  laws  In 
such  a  manner  as  to  emphasize  their  similarity.  Basically,  the  fora  of  the 
equations  Is  fixed  by  the  choice  of  the  coordinate  system  in  which  the  integra¬ 
tion  is  performed,  with  the  critical  factor  being  the  movement  (rotation)  of 
that  coordinate  system  relative  to  the  inertial  system.  The  use  of  the  equation 
for  the  derivative  of  a  vector  referred  to  a  rotating  frame  (equation  6)  is  then 
of  fundamental  Importance,  and  must  be  understood  in  order  to  gain  an  intuitive 
knowledge  of  the  dynamics  involved.  Many  other  key  factors  to  the  understanding 
of  the  missile  dynamics  equations  have  also  been  discuss,  d  and  a  thorough 
understanding  of  these  will  greatly  enhance  the  practicing  engineer's  insight. 
This  then  forms  the  basis  of  this  report  which  is  intended  to  be  both  a  primer 
for  new  engineers  and  a  reference  to  experienced,  engineers  involved  in  missile 
flight  simulation. 
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APPENDIX 


The  approach  used  in  Section  III,  part  B  can  be  extended  to  obtain  the 
displacement  of  the  missile  referred  to  the  missile  body  frame.  This  is  not 
usually  done  as  the  displacement  is  easier  to  visualize  when  related  to  an 
inertial  or  earth  fixed  frame,  but  will  be  done  here  in  order  to  complete  the 
discussion. 

Equation  46  can  be  written  in  terms  of  the  displacement  vector  "R“  which 
is  referred  to  the  body  frame  by  using  equation  6  twice.  Thus 


(71) 


(72) 


F*  -  mR  (73)  . 

or 


F'/m  -  (R)r  +  w  x  R  4-  2w  x  (R)r  +  w  x  (w  x  R)  (74) 


The  acceleration  referred  tb  the  missile  body  coordinate  frame  is  then  given 
by: 


(R)r  -  F’/m  -  w  x  R  -  2 w  x  (R)r  -  w  x  (w  x  R) 

This  equation  can  be  integrated  twice  to  obtain  the  displacement  "R”  ri 
to  the  missile  body  coordinate  frame.  It  is  interesting  to  note,  howe' 
that  the  first  integration  of  equation  75  does  not  provide  the  veloclt; 
referred  to  the  missile  body  frame*  This  can  be  better  understood  by 
examining  the  equation  for  the  velocity  derivative  referred  to  the  bod; 
which  is: 


(75) 

ferred 

er; 

frame 


From  equation  6,  V  is  found  to  be: 

V  -  (R)r  +  w  x  R  (77) 

Now  writing  equation  76  in  the  form: 

F'/m  -  (V)r  +  w  x  V  (78) 

and  substituting  equation  77  for  V  gives: 

F*/m  ■  (R)r  +  w  x  R  +  2w  x  (R)r  +  w  x  (w  x  R)  •  (79) 

which  is  identical  to  equation  74.  Thus  equations  74'  and  78  are  equivalent 
and  comparing  terms  in  the  equations  reveals  that: 


(V)r  »  (R)r  +  w  x  R  +  w  x  (R)r  (80) 

and 


w  x  V  -  w  x  (R)r,  ♦  w  x  (w  x  R)  (81) 

From  equation  80,  it  is  apparent  that  the  derivative  of  V  does  not  equal  the 
second  derivative  of  R  when  both  are  referred  to  the  moving  coordinate  frame 
axis  system.  Thus  integrating  equations  75  and  76  will  give  different  results 
for  the  velocity. 
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